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We analyse the behavior of the thermal conductivity, k(H), in the vortex state of a quasi-two-dimensional 
d-wave superconductor when both the heat current and the applied magnetic field are in the basal plane. At low 
temperature the effect of the field is accounted for in a semiclassical approximation, via a Doppler shift in the 
spectrum of the nodal quasiparticles. In that regime k(H) exhibits twofold oscillations as a function of the angle 
between the direction of the field in the plane and the direction of the heat current, in agreement with experiment. 



Experiments which show that the supercon- 
ducting order parameter in the high-T c cuprates 
is d-wave are often sensitive to the surface effects, 
and there is still interest in the bulk probes of the 
symmetry of the gap. One piece of evidence for 
the linear nodes in the bulk comes from the veri- 
fication of the universal low temperature limit of 
the in-plane thermal conductivity, |Q another is 
based on the observed non-linear dependence of 
the electronic specific heat on the applied mag- 
netic field, H, in the mixed state. j| 

The latter result is based on the observation 
that the properties of a d-wave superconductor 
in the dilute vortex regime (H <C H C 2) are de- 
termined by the near-nodal quasiparticles in the 
bulk.|| In a semiclassical treatment, the energy 
of quasiparticle with the momentum k„, where 
n labels a node, at point r is shifted by Slu = 
v s (r) • k„, where v s is the velocity field associ- 
ated with the supercurrents. In the regions of 
the Fermi surface near the nodes where this shift 
exceeds the local gap, there exist unpaired quasi- 
particles. Since the typical supermomentum is 
Ti/R where R 2 w ^q/ttH, and since the number 
of vortices n v oc H, the spatially averaged density 
of states in a pure c£-wave material varies as \/~H. 
If the positions of the impurities and the vortices 
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are uncorrelated, the argument can be general- 
ized to include the impurity scattering, which de- 
pends on the local density of states. [Q This simple 
approach has worked remarkably well in describ- 
ing the low-temperature thermal and transport 
properties of the vortex state with the field per- 
pendicular to the Cu02 layers. In particular, the 
increase of the T = limit of the in-plane ther- 
mal conductivity with H is well described by the 
semiclassical theory. || 

In relatively three-dimensional cuprates, such 
as YBa2Cu307-5, the same approach applies 
when the field is in the basal plane. In that case 
the magnitude of the Doppler shift depends on 
the relative orientation of the field with respect 
to the nodes, and the density of states in the pure 
limit exhibits fourfold oscillations as a function of 
the angle between the direction of the field and 
the crystalline axes.Q Experimental verification 
of this prediction has not been possible so far and 
is hindered by the extrinsic contributions to the 
specific heat[[z| and the orthorhombicity of the 
material. Q On the other hand, the angular de- 
pendence of the thermal conductivity in the vor- 
tex state has already been observed, J|,[l(| and 
here we analyse it in the semiclassical framework. 

As in Refs.^D we assume a cylindrical Fermi 
surface and a <i-wave gap, and approximate the 
superflow by the velocity field around a sin- 
gle vortex. The field and the thermal gradi- 
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ent are applied at angle a and s to the a axis 
respectively. We consider the regime T, 7 <C 
E H « A , where E H « (« / /2)( 7 r$ A Qb / J ffA c ) 1 / 2 
is the average Doppler shift, and 7 is the low- 
energy scattering rate. Following Refs.|||| we 
obtain that the local change in k is 5k(p)/koq = 
E(a) sin 2 13/ p 2 , where «oo/T — 7rA v 2 /6A is 
the universal thermal conductivity Nq is the 
normal state density of states, A is the gap 
amplitude, Vf is the Fermi velocity, (3 is the 
winding angle of the vortex, p is the distance 
from the center of the vortex normalized to R, 
E(a) = (7r£ ,2 / /8rA )max(sin 2 a,cos 2 a), and L 
is the bare scattering rate. The local fc(r) has to 
be spatially averaged to obtain the field depen- 
dence. When the heat gradient, VT, is parallel 
to the field k\\(H) = (k(t)), where the brackets 
denote the average over a unit cell of the vortex 
lattice. For other relative orientations of VT and 
H the averaging procedure is not clear; it was ar- 
gued[[| that n±(H) — is appropri- 

ate for VT _L H. We therefore take here a simple 
approach of averaging independently the compo- 
nents of the heat current along and normal to the 
vortex, and expect that this procedure gives at 
least qualitatively correct results. Then the lon- 
gitudinal and the Hall thermal conductivity are 
given by hi — k\\ cos 2 (a — e) + k± sin 2 (a — e) and 
K2 = (1/2)|(k^ — K||)sin2(a — e)| respectively, 
with K|| = Koo[l + E (a) Ih(Aq/ Eh)] and 
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1/I + E(a) -E(a) sinh" 



■(1) 



The result for k\ is in agreement with the twofold 
pattern of Ref.|o| at T = 0.8K with e = tt/2, 
see Fig.l. The minima of k.\{q) correspond to 
increased scattering by the vortices when the heat 
current is normal to the field. 

In the regime where T > E H Yu et al. §] have 
measured n± — {k\ ± K2)/V2 with e = 0, found 
a twofold pattern, and explained it as a conse- 
quence of Andreev reflection of quasiparticles in 
the presence of supercurrents; we note that the 
angular dependence of k± is very similar to what 
we obtain at Eh T. In the same regime Aubin 
et al. |l0| observed a fourfold symmetry of K\{a), 
consistent with the picture in which the scatter- 
ing of quasiparticle by the vortices becomes more 
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Figure 1. Angular dependence of n\ for Eh = 
0.05A and T = 0.001A . 



important at high T|Tl]] and this scattering has 
the same symmetry as the gap.|^Jll|] The work 
to explain in detail the high T dependence is in 
progress, and will be reported elsewhere. 
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